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ABSTRACT 

We present the physical properties of [Oiii] emission line galaxies at z > 3 as the tracers of active 
galaxies at IGyr before the peak epoch at z 2. We have performed deep narrow-band imaging 
surveys in the Suhaiu/XMM-Newton Deep Survey Field with MOIRCS on the Subaru Telescope and 
have constructed coherent samples of 34 [Oiii] emitters at z = 3.2 and 3.6, as well as 107 Ha emitters 
at z = 2.2 and 2.5. We investigate their basic physical quantities, such as stellar masses, star formation 
rates (SFRs), and sizes using the publicly available multi-wavelength data and high resolution images 
by the Hubble Space Telescope. The stellar masses and SFRs show a clear correlation known as the 
“main sequence” of star-forming galaxies. It is found that the location of the main sequence of the 
[Oiii] emitters at z = 3.2 and 3.6 is almost identical to that of the Ha emitters at z = 2.2 and 
2.5. Also, we investigate their mass-size relation and find that the relation does not change between 
the two epochs. When we assume that the star-forming galaxies at z = 3.2 grow simply along the 
same main sequence down to z = 2.2, galaxies with M* = IO^-IO^^Mq increase their stellar masses 
significantly by a factor of 10-2. They climb up the main sequence, and their star formation rates also 
increase a lot as their stellar masses grow. This indicates that star formation activities of galaxies are 
accelerated from z > 3 towards the peak epoch of galaxy formation at z ~ 2. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: high-redshift 


1. INTRODUCTION 

The activities of star formation in galaxies and those 
of active galactic nuclei (AGN) are very high at z ^ 
1 — 3, corresponding to about 8-10 billion y ears ago 
(e.g. iHookins fc Beacomll2006l iFan et all 120041 1. Physi¬ 
cal states of galaxies are expected to change dramatically 
during this epoch, and it is critical to investigate galaxy 
properties in detail at this epoch in order to under¬ 
stand physical processes of galaxy formation. With the 
near-infrared (NIR) photometric and spectroscopic ob¬ 
servations with ground-based and space telescopes, stud¬ 
ies of high-z galaxies have advanced substantially in re¬ 
cent years, and our knowledge of physical states of those 
galaxies have been expanded significantly by many pre¬ 
vious studies (e.g. lErb et al.l[^06al: Ivan Dokkum et al.l 
20081: iKr iek et al] I2009bt iFbrster Schreiber et al.l 1200ft 
Wuvts et al.ll20lih . 

An example is the discovery of the relationship be¬ 
tween stellar mass and star formation rate (SFR) of 
star-forming galaxies both at low and high redshifts. 
Many previous studies have revealed that these two 
quantities show a tight correlation, ca lled the “main se¬ 
quence” of sty-forming galaxies fe.g-. lElbaz et al.l[2007t 
iNoeske et al.l l2007t iDaddi et al.1 l2007t IWhitaker et al.l 
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120121 : iKashino et al.1 l2013f l. Using galaxy samples 
across a wide redshift range, the evolution of this rela¬ 
tion has been invest i gated (e.g. IWhitaker et al.l I2012t 
iKovama et al.l 120131 : iTasca et al.l 1201^0 They have 
shown that the M*-SFR relation does evolve strongly 
with redshift at least up to z ^ 2.5, in the sense that SFR 
increases monotonically with redshift at a given stellar 
mass. 

Another example is the discovery of compact, mas- 
sive, and quiescent galaxies at z ^ 2 (“red nuggets”; 
iDaddi et al.ll2005HDamianov et al.l2009ll . and their likely 
progenitors, namely, compact star-forming galaxies at 
z > 2. The compact star-forming galaxies are thought to 
be formed by gas-rich processes such as mergers or disk 
instabilities which invoke starbursts in the central com¬ 
pact regions. They would then become compact, qui¬ 
escent gal axies when their s tar formation activities are 
quenched (|Barro et al.l 120131) . We can search for such 
compact star-forming galaxies based on the stellar mass 
versus size diagram. The redshift evolution of the mass- 
size relation is crucial for identifying the evolutionary 
paths from star-forming galaxies to quiescent galaxies, 
and it has been investigated up to z ^ 3. It is found that 
the sizes of early-type galaxies depend more strongly on 
their stellar masses as compared to late-type galaxies and 
that the average size evolution at a fixed stellar mass of 
late-type galaxies is very slow, while the average size of 
the early- type galaxies increases rapidly with the co smic 
time (e.g. iShen et ^120031 : Ivan der Wei et al.ll2014ll . 

The studies of physical properties of galaxies are 
now being expanded to more distant Universe beyond 
the highest peak of the star formation activities at 
z ^ 2 — 3. The epoch of z ^ 3 — 3.7, correspond¬ 
ing to about 1-2 Gyr before the peak epoch, is es- 
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pecially crucial to reveal how galaxy formation activi¬ 
ties are activated towards its peak. At z > 2.5, the 
ultraviolet (UV) light is often used to construct star¬ 
forming galaxy samples. Using the UV-selected galax¬ 
ies, such as the Lyman Break Galaxies (LBGs), the red- 
shift evolution o f star formation activities have been in¬ 
vestigated (e.g. iStark et ^ 120091: iGonzalez et ahll^OlOl : 
Reddy et aj.l 120121 iStark et al.ll2013t iTasca et al.l 1201411 . 

Stark et al.l (1201311 investigate the evolution of the spe¬ 
cific star formation rate (sSFR = SFR/M*) for the spec¬ 
troscopically confirmed LBG sample at 3.8 < z < 5. 
Their result indicates that the sSFR at the fixed stellar 
mass increases by a factor of 5 from z ^ 2 to z ^ 7. 
Based on the i-band s elected and s pectr oscopically con- 
hrmed galaxy sample, iTasca et al.l (1201411 investigate the 
evolution of the M*-SFR and M*-sSFR relations up to 
z ^ 5. Their conclusion is that the sSFR increases very 
slowly from z ~ 3 up to z ~ 5. These two galaxy sam¬ 
ples have some small differences in the sample selection 
(see the papers for details). Moreover, both of the sam¬ 
ples are likely to be biased to less dusty star-forming 
galaxies since the samples are selected at the rest-frame 
far UV. Up to z 2.5, the Ha emission line is a good 
tracer of star formation activity, due to its lower sen¬ 
sitivity to the dust extinction as compared to the UV 
light. However, since the Ha line comes to longer wave¬ 
length than A-band at z > 3, we are no longer able to 
use the Ha line as the tracer of star-forming galaxies at 
z > 3 in the observations with ground-based telescopes. 
Recent NIR spectroscopic observations have revealed 
that the high-z star-f orming galaxies show strong [Oiii] 
emission lines (e.g. [Masters et al.l 12014 iHolden et al.l 

12014 [Seidel et all 2014 iShimakawa et al.l[2014ll . Such a 

strong [Oiii] emission line indicates extreme interstellar 
me dium (ISM) conditions of high-z star-forming galax¬ 
ies. iNakaiima fc Ouchil (j2014f ) have shown that the lower 
metallicities and higher ionization parameters contribute 
to their extreme ISM conditions. Therefore, we argue 
that the [Oiii] emission line is one of the best tracers of 
star-forming galaxies at z > 3. We should note however 
that the [Oiii] emission line originates from ionized re¬ 
gions not only by hot young massive stars but also by 
AGNs. 

We have conducted a systematic narrow-band imag¬ 
ing survey with Subaru Prime Focus Gamera (Suprime- 
Cam) and Multi-Object InfraRed Camera and Spectro¬ 
graph (MOIRCS) on the Subaru telescope. The project 
is called “MAHALO-Subaru ” (MApping HAlph a and 
Lines of Oxygen with Subaru; iKodama et al.l[201^ . The 
imaging observations with narrow-band filters allow us 
to obtain emission line galaxies in particular redshift 
slices. By targeting high density regions, such as clus¬ 
ters or protoclusters, as well as the lower density blank 
helds, we have constructed a coherent NB-selected galaxy 
sample across vari ous environments and cosmic times. 
iTadaki et al.l (120131 1 carried out NB imaging observations 
with MOIRCS at the Subaru/AMM-A ewton Deep sur¬ 
vey Field ("SXDF: fFunisawa et al.l[2f)?)^ . and have con¬ 
structed the Ha emitter sample at z=2.2 and 2.5. They 
found that the Ha emitters constitute the star-forming 
main sequence and that galaxies with high SFRs located 
above the main sequence tend to be dustier. They also 
investigated the structural properties of the Ha emit¬ 
ters using the Hubble Space Telescope (HST) images 


(|Tadaki et al.ll20l4l . They found two classes of objects; 
star-forming galaxies with extended disks and massive 
and compact star-forming galaxies. The latter class of 
objects are likely to evolve to the red nuggets seen at sim¬ 
ilar or lower redshifts after quenching their star formation 
activities and eventually to massive elliptical galaxies to- 
day after significant size growth by minor mergers, (e.g. 
Ivan Dokkum et ^ 1201011 In order to make such evolu¬ 
tionary links further, it is important to explore those 
possible progenitors (massive and compact star-forming 
galaxies) at higher redshifts, z > 3, in a systematic way. 

In this study, we construct a NB-selected [Oiii] emitter 
sample at z > 3 based on the NB imaging survey data 
at SXDF taken through the MAHALO-Subaru project 
(|Tadaki et al.ll2013ll . We investigate their physical prop¬ 
erties to reveal galaxy formation at the epoch shortly be¬ 
fore the highest peak of activities. In this field, a range 
of multi-wavelength photometric data and high resolu¬ 
tion images by HST are available. Combining those data 
with our NB imaging data, we investigate stellar masses, 
SFRs, and sizes of the [Oiii] emitters. We will then 
compare their properties with those of the Ha emitters 
at lower redshifts (2.2 < z < 2.5) in the same held to 
track the evolution of star-forming galaxies between the 
two epochs. 

This paper is organized as follows: In Section 2, we 
briehy introduce our NB imaging observations and the 
publicly available observational data in SXDF and ex¬ 
plain the selection method of emission line galaxies. In 
Section 3, we present the measurements of basic physical 
quantities of the [Oiii] emitters, such as stellar masses, 
SFRs, and sizes, and show the relations between these 
physical quantities. We also compare our galaxy sam¬ 
ple with the Ha emitters at lower redshifts in the same 
Held. In Section 4, we discuss any possible selection bias 
arising from our use of the [Oiii] emission line as the 
tracer of star-forming galaxies at z > 3. We also discuss 
the mass growth and star formation histories of galaxies 
from z = 3.2 to z = 2.2 with a simple model. Finally, 
we summarize our study in Section 5. We assume the 
cosmological parameters of Dm = 0.27, Da = 0.73, and 
Hq = 71 [km s“^ Mpc“^]. Through this paper, unless 
otherwise noted, al l the magnitudes ar e given in the AB 
magnitude system (lOke fc GunnllI983ll , and the Salpeter 
initial mass function ( SalDeterllI955ll is adopted for the 
estimations of stellar masses and SFRs. 


2. OBSERVATIONAL DATA AND SAMPLE SELECTION 
2.1. Narrow-band Imaging Surveys 

We summarize below our NB imaging observations per¬ 
formed at SXDF, and further detai ls of the observation s 
and data reduction are described in ITadaki et al.l (j20I3ll . 

The NB imaging at SXDF was carried out with 
the NIR camera and spectrograph called MOIRCS 
([Suzuki et al.ll200§l on the Subaru telescope. MOIRCS 
is equipped with two HAWAH-2 detectors (2048 x 2048 
pixels). The pixel scale is 0.II7 arcsec/pixel and the 
Held-of-view (FoV) is 4' x 7'. Two NB hlters were used, 
namely, NB209 (Ac = 2.093 ^im, FWHM = 0.026 urn) 
and NB23I5 (Ac = 2.317 /im, FWHM = 0.026 /rm). 
The NB209 and NB23I5 Hlters can probe Ha emission 
lines at z = 2.I9I ± 0.019 and z = 2.525 ± 0.021, and 
also [Oiii]A5007 emission lines at z = 3.174 ± 0.025 and 
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Fig. 1.— Transmission curves of the two narrow-band filters on 
MOIRCS; NB209 and NB2315 (solid lines), and the .R-band fil¬ 
ter on UKIRT/WFCAM (dotted line). X-band data is used to 
estimate the stellar continuum (Section 2.3). 

z = 3.623 ± 0.027, respectively. Figure [T] shows the 
transmission curves of these two NB filters. The data 
were obtained over several observing runs from October 
to November in 2010 and in September 2011. The to¬ 
tal observed areas were 91 arcmin^ and 93 arcmin^ for 
NB209 and NB2315, respectively. The survey areas are 
slightly different between the two NB filters, but both 
are overlapped with the CANDELS-UDS field observed 
by H ST/AGS and WFC3 (see Figure 2 of iTadaki et al.l 
l2013f) . The exposure times were 140-186 minutes, and 
the seeing sizes were 0.5" — 0.7" (FWHM). The 5cr lim¬ 
iting magnitudes with 1.6" diameter aperture were 23.6 
mag and 22.88 mag in NB209 and NB2315, respectively. 
The observed data were reduce d with the MOIRCS imag¬ 
ing pipeline software iMCSRED^ dTanaka et al.l[2011[l . The 
point spread functions (PSF) were smoothed to 0.7" 
when combining all the images. 

2.2. Public Data 

In our survey field, SXDF, the multi-wavelength data 
from UV to mid-infrared (MIR) are all available. We 
use the public phot ometric catalog provi ded at the 
Rainbow Database^ (iGalametz et al.l l2013f) . It con¬ 
tains u-band data from CFHT/Megacam (O. Almaini 
et al. in preparation), B-, V-, Rc-, i '-, and z'-band 
data from Subaru/Sprime-Cam ISXDS: [Furusawa et al.l 
EOOl), Y- and A:,-band d ata from VLT/HAWK-1 
(HUGs; iFontana et al.l I2014H . J-, H- , and jF-band 
Ata from UKIRT/WFCAM fUKlDSS: iLawrence et al.l 
I2007I1 . 3.6/im, 4.5/rm, 5.8/rm, and 8.0/im data from 
Spitzer/IRAC, and 24 fj,m data from Spitzer/Ml PS 
(SpUDS; PI: J. Dunlop and SEDS: lAshbv et aT1l2013ll . 

Our survey are as are also mostly covered by the 
ffST/ CANDELS (|Grogin et al.l 1201 it iKoekemoer et al.l 
1201 If) , and the photometric catalog and the high reso¬ 
lution images of Vpeoew-, 7_F8i4W-band from ACS, and 
Jfi 25 W-, 77Fi60iv-band from WFC3 are all publicly 
available. All the publicly available multi-wavelength 
data used in this study are summarized in Table [TJ 


TABLE 1 
SXDF Data Set 


Instrument 

Filter 

5(7 limiting magnitude 

CEHT/MegaCam 

u 

27.68 

Subaru/Suprime-Cam 

B 

28.38 


V 

28.01 


Rc 

27.78 


i' 

27.69 


z' 

26.67 

HST/ACS 

F606VF 

28.49 


F814VF 

28.53 

HST/WFC3 

Fl2bW 

27.35 


Fimw 

27.45 

VLT/HAWK-I 

Y 

26.73 


Ks 

25.92 

UKIRT/WECAM 

J 

25.63 


H 

24.76 


K 

25.39 

Spitzer/IK AC 

3.6pm 

24.72 


4.5pm 

24.61 


5.8pm 

22.30 


8.0pm 

22.26 

Spitzer/M.IPS 

24/rm 

30-60 pjy 


K- bands) images taken with WFCAM, using th e pub¬ 
lic software SExtractor (jBertin fc Arnout^ll996[ ). The 
pixel scales and PSF sizes of the BB images are matched 
to those of the NB images. We perform aperture pho¬ 
tometries on the NB and BB images with a 1.6" diam¬ 
eter aperture. Source extraction and photometries are 
carried out with the double image mode of SExtractor. 
Those aperture photometry data are used to select line 
emitters based on the color-magnitude diagrams (Figure 
HI) and to measure NB fluxes. On the other hand, the 
tem plate-fitting photomet ry data from the public catalog 
(see IGalametz et ^120131 for more details) are used for 
the color-color selections (Figure [3]), the spectral energy 
distribution (SED) fitting, and the SFR measurements 
from UV luminosities (in i?c-band). 

The objects that have large excesses in NB fluxes as 
compared to BB fluxes are selected as the NB emitters. 
Since the effective wavelengths of the NB filters and the 
AT-band are slightly different, we estimate BB fluxes at 
the exact effective wavelengths of the NB209/NB2315 
filters by i nterpolating fluxes between H- and AT-bands 
as follows (ITadaki et a.l.ll201^ : 

HK{X = 2.09 fim) = 0.8 K + 0.2 H - 0.015, (1) 

HK{\ = 2.315 Aim) = 1.2 AT - 0.2 H + 0.011. (2) 

We select NB excess sources using FfAT—NB209 or HK— 
NB2315 color-magnitude diagrams (Figure [2]). We use a 
parameter E, which determines the signifi cance of a NB 
exces s relative to an Icr photometric error (jBimker et al.l 
Il995f) . The relation between E and the color of rriBB — 
tunb is obtained from 


2.3. Selection of [Oiii] emitters at z > 3 

First of all, we extract sources from the NB images 
taken with MOIRCS and the broad-band (BB; H- and 


ruBB — ruNB = —2.5 logio 


, _ ^a/'^BB ii'NB 
/nb 


(3) 


1 http://www.naoj.org/stafr/ichi/MCSRED/mcsred.html where /nb is a NB flux density, crlg and cr^B ^he 

^ https://arcoiris.ucolick.org/Rainbow_navigator_public/ sky noises in the BB and NB images measured within an 
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aperture, respectively (|Tadaki et alJl201^ . We set the 
criterion of E = 3 to sample secure emitters, and it is 
represented as the solid curve in Figure [2j We also set 
the criteria that the NB magnitude excess with respect 
to the HK magnitude, HK — NB, is larger than 0.4 
mag, and that the NB magnitude is brighter than the 
5(7 limiting magnitude (the horizontal red line and the 
vertical dashed line in Figure^ respectively). With these 
selection criteria, we obtain 101 and 58 NB excess sources 
in NB209 and NB2315, respectively. 

The NB209(NB2315) filter captures different emission 
lines from galaxies at different redshifts, such as Ha 
emission lines at z = 2.19(2.53) and [Oiii] emission lines 
at z = 3.17(3.62). In order to construct a sample of 
[Oiii] emitters at z > 3, we need to separate [Oiii] emit¬ 
ters from Ha emitters at z < 3. The conspicuous spec¬ 
tral features such as the Balmer and/or 4000A breaks 
can be used for this purpose. Since the break feature 
is redshifted to a different wavelength in the observed 
frame according to the redshift of a galaxy, if we choose 
an appropriate combination of passbands that can neatly 
straddle the spectral break feature, we are able to disen¬ 
tangle those various possible s olutions for di f ferent line 
emitters at different redshifts. iTadaki et al.l (|2013f) use 
i' — J versus J — K diagram for NB209 emitters and 
i' — H versus H — K diagram for NB2315 emitters. They 
have shown that [Oiii] emitters at z > 3 and Ha emit¬ 
ters at z < 3 can be well separated by the dividing lines 
as shown in Figure [31 We follow their selection method 
and here the contribution from the emission line is sub¬ 
tracted from the if-band magnitude. Our NB emitter 
sample may also contain some Pa emitters at low red¬ 
shifts (z =0.1-0.2). However, since the Pa emitters 
should appear in the same regions as the Ha emitters on 
these diagrams, they should not be major contaminants 
for our [Oiii] emitters. We have thus finally obtained 
strong candidates for [Oiii] emitters; 27 at z = 3.17 and 
7 at z = 3.62. 

Figure m shows the relation between stellar masses and 
dust-extinction-uncorrected SFRs measured from [Oiii] 
line luminosities for the [Oiii] emitters in order to ver¬ 
ify our sample selection. We will explain our method 
of measuring stellar masses and SFRs in the following 
sections. Our criteria of determining NB flux excesses, 
namely E > 3 and HK — NB209 > 0.4 [mag] as shown in 
FigurejU correspond to the limits of SFR > 4.5 [M0yr“^] 
and EWrest > 30 [A], respectively. We draw a line corre¬ 
sponding to the EW cut on the M,-SFR diagram by es¬ 
tablishing a relation between stellar mass and SFR along 
the threshold of HK — NB209 = 0.4. As proxies of the 
H — K and J — K colors along this boundary line, we 
use the averaged colors of the 3 objects which are located 
nearest to the boundary of HK — NB209 = 0.4. Then, 
for a given NB magnitude, we assign a if-band magni¬ 
tude using HK — NB209 = 0.4 and the above H — K 
color. A stellar mass is then estimated from the if-band 
magnitude and the above J—K color. J—K color is used 
to correct for the mass-to-light ra tio based on the stel¬ 
lar p opulation synthesis model of iKodama et al.l (jl998l 
119991 1. Also a dust-extinction-uncorrected SF^qiii] is 
calculated from each NB magnitude. Figure shows 
that our EW cut is located well below the actual ob¬ 
served data points, and our selected galaxy sample is not 


biased to any particular galaxies on the main sequence 
diagram. 

3. PHYSICAL PROPERTIES OE [Oiii] EMITTERS 
3.1. AGN Contribution 

Not only hot young massive stars in star-forming re¬ 
gions but also AGN activities can contribute to the [Oiii] 
line intensity. Therefore it is important for us to verify 
the presence of AGN candidates in our sample. 

We use a phot ometric redshift code EAZY 
(jBrammer et al.l 1200811 to obtain the rest-frame [/-, 
V-, and J-band magnitudes for our sample. Note that 
derived photometric redshifts are mostly consistent with 
the expected redshifts, z = 3.2 and 3.6. The rest-frame 
U — V and V — J colors allow us to distinguish between 
two galaxy populations, namely, old quiescent galaxies 
and young, dusty star-forming galaxies, by capturing 
the Balmer/4000A breaks between U- and F-bands (e.g. 
Wuvts et al.ll2007l : I Williams et al.ll200l : iWhitaker et al.l 
201111 . Figure [5l shows the rest-frame UV J color-color 
diagram for our [Oiii] emitters. We find that one object 
is marginally classified as a quiescent galaxy, indicating 
that its [Oiii] emission is likely to be dominated by the 
AGN activity rather than the star formation. However, 
we cannot discriminate between the contribution from 
AGNs and that from star-forming regions for all the 
other emitters classified as star-forming galaxies. 

For further investigatio n, we inspect the X-ray image 
by the XMM — Newton (lUeda et al.ll2008il . None of our 
[Oiii] emitters is detected in X-ray and thus our sample 
does not seem to contain any bright unobscured AGNs. 
We also look into the Spitzer/MAPS 24^m catalog and 
find that three objects are detected with MIPS. A frac¬ 
tion of them might be obscured AGNs with warm dust 
components which emit strong IR emissions. 

Spectroscopic observations are necessary to confirm 
the presence of AGNs and we do not exclude these ob¬ 
jects in the following analyses. 

3.2. SED fitting 

We perform the SED f itting for our [Oiii] emitters us¬ 
ing a public code FAST (|Kriek et al.l 12009^ . We use 18 
bands, u, B, V, Rc, i', z', F606W, F8141F, F125VF, 
F160W, Y, J, H, Kg, 3.6^m, 4.5^m, 5.8^m, and 8.0^m. 
For the NB209-selected [Oiii] emitters, emission line 
fluxes are subtracted from the ATg-band fluxes before the 
SED fitting is performed, while no correction is required 
for the NB2315-selected ones as the NB2315 has little 
overlap with the ATg-band in wavelength. The redshifts 
of the NB209- and NB2315-selected [Oiii] emitters are 
fixed to z = 3.17 and z = 3.62, respectively, for the SED 
fitting. We use the s tellar population synthesis model of 
Bruzu al &: Gharlotl ([200311 . the Salpeter IMF (ISalpetei 

S and the dust attenuation law of iCalzetti et al. 

. We assume the exponentially declining star for¬ 
mation history in the form of SFR exp(—t/r) with 
log(r/yr) = 7.0-10.0 in steps of 0.5, and the solar metal- 
licity. The output physical quantities from the FAST code 
are star formation timescale r, age, dust extinction Ay, 
stellar mass, SFR, specific SFR and age/r ratio. 

The stellar masses and dust extinctions (Ay) used in 
the following analyses are all estimated by the SED fit¬ 
ting. 
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Fig. 2.— The color-magnitude diagrams for NB209 (left) and NB2315 (right). Gray dots represent all the sources detected in the NB 
images, and blue filled circles are the sources selected as the NB emitters with our selection criteria. Red open circles show the [Olll] 
emitters at z = 3.17 (left) and z = 3.62 (right) selected based on the color-color diagrams (FigurelSjl. The solid curves correspond to ±3 (t 
photometric errors and the horizontal red lines correspond to HK — NB = 0.4 [mag], which are used to select the NB emitters. The vertical 
dashed lines correspond to the 5a limiting magnitudes of the NB images. 




i — J [mag] i - H [mag] 

Fig. 3.— The color-color diagrams for NB209 (left) and NB2315 (right). Blue filled circles represent the NB emitters selected based on 
the color-magnitude diagrams (Figuref^, and red circles show the [Olll] emitters at z = 3.17 (left) and z = 3.62 (right) which satisfy the 
criteria shown by the dividing solid lines. The arrow represents the reddening vector of Av = 1 mag. 


3.3. Star Formation Rates 

We estimate SFRs of the [Oiii] emitters with the two 
different indicators, namely, the UV continuum luminosi¬ 
ties (tracing hot young stars) and the [Oiii] emission 
line intensities (tracing star-forming HII regions). In 
the former case, w e adopt the following equation from 
iMadau et al.l (|1998l) : 


SFR (M —_ 

(1 + z)8 X 10^’^ (erg cm^^s^^Hz”!) 

L(1600A) 

8 X 10^”^ (erg s“^Hz“^)’ 

where is the luminosity distance and /i, is the flux 
density derived from the i?c-band magnitude (Ac = 
6498. lA). The dust extinction at 1600A is estimated from 
the SED-based Ay value and the extinc tion curve for 
starburst galaxies of iCalzetti et al.l (l2000[ ): 

EiB - VUellar = A'y/R'y, (5) 


Al'(A) = k'{X)EiB - E)stenar, (6) 


fc'(A) = 2.659(-2.156 -h 1.509/A - 0.198/A2 -h O.Oll/A^) 
-|-i?Y, at 0.12^m < A < 0.63^m. (7) 

E{B — E)stenar indicates the amount of reddening in the 
stellar continuum, a nd i?(, is 4.05 for starburst galaxies 
(|Calzetti et al.l[2000l) . The intrinsic flux density fi{X) is 
then obtained as 

/AA) = /o(A)10°-4-^'(^\ (8) 

and the dust-extinction-corrected SFRs (SFRuv) are de- 
rive d using Eq([4l). 

In lMaschiettoet "m (|2008[) . they derive a SFR from an 
[Oiii] emission line strength by assuming the [Oiii]/Ha 
ratio of ^ 2.4 which is the maximum valu e for local 
star-forming galaxies (|Moustakas et al.ll2006D . Consid¬ 
ering the fact that high-z star-forming galaxies show 
very high [O iii]/H;9 ratio due to the high exc i tation 
states fe.g. Ma sters et al.1 12014 iHolden et al.l 12014 
iSteidel et al.ll20l4lShimakawa et al.ll2014D . this assumo- 
tion seems to be reasonable for our sample, although this 
ratio has a large dispe rsion among individual galaxies 
([Moustakas et al.1l2006D . We adopt this maximum ratio 
to the relation be tween the SFR and Ha luminosity of 
iKennicufrB (|1998bD : 


and 















































6 


Suzuki et al. 



M, [Mg] 

Fig. 4.— The stellar mass versus SFR diagram for the [Olll] 
emitters before dust correction is applied to discuss any possible 
selection bias. Filled and open circles represent the [Oiii] emitters 
at z = 3.2 (NB209) and 2 = 3.6 (NB2315), respectively. SFRs 
are derived from [Oiii] line luminosities without dust correction 
here on purpose. The horizontal solid line and the dot-dashed line 
correspond to our selection criteria of the NB emitters, namely, E = 
3 and HK — NB = 0.4, respectively. These numbers correspond 
to our detection limits of SFR^qiii] = 4.5 [M 0 yr~^] and EWrest = 
30 [A], respectively. 



V — J [mag] 

Fig. 5.— The rest-frame UV J diagram for our [Oiii] emitters 
(filled circles for z = 3.2 and open circles for 2 = 3.6). Star marks 
indicate MIPS detected sources. The arrow shows a reddening 
vector of E{B — y) = 0.2. The typical errors of the U — V and 
V — J colors are shown at the bottom-right corner of the figure. 
We make 1000 pseudo photometric catalogs for each object by ran¬ 
domly generating fiuxes in each band from a Gaussian distribution 
with the standard deviation determined by each flux error. These 
generated SEDs are fitted by EAZY and the dispersions of the two 
colors are estimated for each object. We show those dispersions as 
the typical error of each color. 


SFRho (Mgyr-i) = 7.9 x (9) 

erg 

The lower limit of SFR^qiii] is thus obtained by; 

SFR[oiii] (Meyr-i) > 0.33 x (lo) 

J gj.g g i 

The luminosity of the [Oiil] emission line, T[oiii], is ob¬ 
tained by measuring the [Oiil] line flux from the NB and 



Fig. 6. — The comparison between the dust extinction corrected 
SFR^OIII] SFRuv (filled circles for 2: = 3.2 and open circles 
for 2 : = 3.6). Arrows represent that SFR^qiii] is shown as the 
lower limit. Star marks show MIPS detected objects, and a square 
indicates an object classified as a quiescent galaxy based on the 
UV J diagram (Figure 0 . The dashed line represents the case 
where two measurements are identical. 


BB flux densities. The NB and BB flux densities are de¬ 
fined as 


/nB = fc + Fjine/ANB, 

(11) 

/bB = /c + Fjine/ABB, 

(12) 

where fc is a continuum flux density, Fjine is a line flux 
intensitv. and Anr and Arr are FWHMs of the NB 
and BB filters (|Tadaki et al.[l20I3fl. The continuum flux 
density, the line flux intensity, and the equivalent width 
(EW) in the rest-frame are given by the following equa¬ 
tions, respectively; 

r _ /bb —/nb(Anb/Abb) 

I — Anb/Abb 

(13) 

. /nB - /bb 

Fiine = ANB * , 

i — zanb/Abb 

(14) 

EWrest =^““(1 + ^)-^ 

Jc 

(15) 


The line flux Fune is converted to the [Oiii] luminos¬ 
ity with T[oiii] = 47rn|^Fiine. The dust extinction at 
5 OO 7 A is estimated in the same manner as used for the 
dust extinction at I 6 OOA based on the SED fitting. We 
assume that there is no extra extinction for the neb¬ 
ula emissions compared to the stellar extinction, i.e. 

E{B — R)stellar = E(i3 — R)nebula- 

In Figureini we compare SFRs derived from the two dif¬ 
ferent indicators. The ratios of SFR[oiii]/SFRuv range 
from 0.25 to 3 for most of the objects. On the other 
hand, the object classified as a quiescent galaxy on the 
UVJ diagram (Figure [S|) shows a slightly higher ratio of 
SFR[ohi]/SFR uv ~ 5-3. It suggests that this object has 
an extra contribution from an AGN to its [Oiil] emission 
as expected in Section 3.1. 
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3 . 4 . M^-SFR Relation 

We investigate the relation between stellar masses and 
SFRs (the “main sequence” of star-forming galaxies) 
for the [Oiil] emitters at z > 3. The dust-extinction- 
corrected SFRs (SFRuv) of most of the [Oiil] emitters at 
SXDF range from a few to ^ 30 MQyr“^. Figure 0 shows 
a M*-SFR relation for the [Oiil] emitters at z=3.2 and 
3.6, together w ith the Ha emitter s at z=2.2 and 2.5 in 
the same field (iTadaki et al.ll201^ . We re-estimate stel¬ 
lar masses and dust-extinction-corrected SFRs (SFRuv) 
of the Ha emitters in the same manner as in this study 
by using the FAST code (Section 3.2). 

We find that stellar masses and SFRs of the [Oiil] 


emitters show 
studies of star 
shift range (e.g 

a clear correlation as seen in other 
-forming galaxies across a wide red- 

. lElbaz et al.l 12007 

: IDaddi et all 120071: 

Whitaker et al.1 

20I2HKovama et all 

20131: IKashino et al.l 

20131: lTa.sca et alJl2ni4l). The normalization of the M*- 


SFR relation for the [Oiil] emitters at z=3.2 and 3.6 
looks almost identical to that of the Ha emitters at 
z=2.2 and 2.5. We confirm that the best-fit line to the 
[Oiil] emitters is consistent with the fit to the Ha emit¬ 
ters within Icr errors in both the slopes and the inter¬ 
cepts. Importantly, however, the distributions of galax¬ 
ies along the main sequence are systematically different 
in the sense that the stellar masses of the [Oiii] emitters 
at z ^ 3.2 and 3.6 are much (nearly by a factor of 10) 
lower than those of the Ha emitters at z ^ 2.2 and 2.5. 
We can simply interpret the difference in two ways; (1) 
the evolution of galaxies from z=3.2, 3.6 to z=2.2, 2.5, 
and (2) the selection bias between [Oiii] and Ha emit¬ 
ters. In Section 4.1, we will refer to the option (2), and in 
Section 4.2, we will assume that the difference between 
stellar mass distributions is only due to the evolution of 
galaxies from z=3.2, 3.6 to z=2.2, 2.5, and discuss how 
their stellar masses and SFRs grow in this time interval. 

3 . 5 . Sizes 

Our NB imaging survey areas are covered by the 
77ST/CANDELS fields, and high spatial resolution im¬ 
ages in the rest-frame UV-optical wavebands (ACS and 
WFC3) are available for the [Oiii] emitters. 

We use the s t ructu ral parameters measured by 
Ivan der Wei et all (|2012f) on the Ffuieow-band selected 
objects in CANDELS. We briefly summarize below their 
methods to obtain the structural parameters. They per¬ 
form a Sersic model fit to the Hpifin w-band selected ob- 
jects using t he softwares; GALA PAGOS (|Barden et al.ll201^ 
and GALFIT (|Peng et al.ll201Clll . The parameters that are 
used for the fit are the total magnitude, half-light radius 
measured along the major axis, Sersic index, axial ratio, 
position angle, and central position. The initial guesses 
of these parameters are given by SExtractor. The 
best-fit GALFIT parameters for all the objects are pub¬ 
licly available in the Rainbow Database (iGalametz et al.l 
l20I3f) . A flag number between 0 and 3 is assigned for 
each object. We reject objects with flag > 2 because the 
fitting result with a Sersic model becomes increasingly 
unreliable. They note that resultant structural parame¬ 
ters have systematic and random uncertainties and that 
such uncertainties depend on the brightness of objects 
and become larger for fainter objects. At Hfiqow = 25, 
systematic and random uncertainties of the half-light ra- 



Fig. 7. — The stellar mass and SFR relation (main sequence of 
star-forming galaxies) for the [Oiii] emitters (filled circles for 2 : = 
3.2 and o pen circles for z = 3 .6) and the Ha emitters at 2=2.2 and 
2.5 from H^daki et all 1120131') (filled triangles). SFRs are derived 
from UV luminosities and corrected for dust extinction. Errors in 
SFRs for the [Oiii] emitters are estimated from Icr photometric 
errors of the Rc-band magnitudes. The solid line represents the 
best-fit line to the Ha emitters; SFRyv = 129 {Mu = 

M*/10^^Mq). Top and right-side histograms show the stellar mass 
and SFR distributions of the [Oiii] emitters at 2=3.2 and 3.6 (red 
hatched histograms) and the Ha emitters at 2=2.2 and 2.5 (blue 
open histograms). 


dius (re in arcsec) are 0.04" and 0.18", respectively, when 
Tg is less than 0.3", while they are —0.09" and 0.33" when 
re is greater than 0.3". At Hpieow = 26, those uncer¬ 
tainties become as large as 0.12"/0.42" (re < 0.3") and 
-0.II"/0.63" (re > 0.3"), respectively^ _ 

Usin g the GALFIT parameters from Ivan der Wei et al.l 
(IMl . we estimate the effective radius [kpc] in the 
rest-frame t/-band for the [Oiii] emitters. We use only 
the bright objects {Hfiqqw < 25) with the flag values 
of 0 or 1. Based on the magnitude cut and the flag 
values, 10 and 9 objects, respectively, are excluded from 
the [Oiii] emitter sample. 

3 . 6 . M^:-Size Relation 

Figure [5] shows the relation between stellar masses and 
sizes for the [Oiil] emitters at z > 3. The Ha emitters 
at z = 2.2 and 2.5 are also shown. Their sizes are esti¬ 
mated in the HST Ji 25 w-band images so that they can 
be directly compared to those of our [Oiii] emitters at 
the same rest-frame wavelength. The size measurements 
are also limited to the bright objects {Jfi 25 W < 25) 
with flag = 0 or I. In Figure (5] the solid and dashed 
lines represent the mass-size relations of late- and early- 
type galaxies at 2.5 < Zphot < 3, re spectively, derived 
from the 3D-HS'T/CANDELS group (|van der Wei et all 
l20I4f l. We find that the size distribution of the [Oiii] 
emitters with respect to the stellar mass is similar to 
that of the Ha emitters, and that they follow the mass- 
size relation of late-typ e ga laxies at z ^ 2.75 from 
Ivan der Wei et al] (1201^ . In Ivan der Wei et al.l (I20I4I1 . 
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Fig. 8 .— The stellar mass and size relation for the [Olll] emitters 
(filled circles: z = 3.2, open circles: z = 3.6) with Hp\^nw < 25 
and the Ha emitters at z=2.2 and 2.5 from iTadaki et ahl 1120131 ) 
(filled triangles) with Jfi 251V < 25. The solid and dashed lines 
represent the mass-siz e relations of 2: 2.75 la te- and early-type 

galaxies, respe ctively llva n der Wei et al .l 120141 ). Note that their 
galaxy sizes in Ivan der Wei et al.l II2014I1 are estimated at a rest- 
frame wavelength of 5000A, slightly longer wavelength than the 
rest-frame 17-band where our galaxy sizes are measured. 

galaxy sizes are estimated at a rest-frame wavelength of 
5 OOOA, slightly longer wavelength than the rest-frame U- 
band where our galaxy sizes are measured. To verify a 
possible effect due to wavelength mismatch, we also ap¬ 
ply their same correction method to our sample, and con¬ 
firm that there is no systematic difference between the 
two measurements at different wavelengths. 

While most of the [Oiil] emitters have sizes consistent 
with the mass-size relation of late-type galaxies at z ^ 
2.75, there is a massive [Oiii] emitter for its size (M* ^ 
3 X and re ~ 1 kpc). Massive and compact 

star-forming galaxies are expected to evolve to massive 
and compact quie scent galaxies when their star formatio n 
is quenched (e.g. iBarro et al.l 120131 : ITadaki et al.l 120141) . 
We confirm the presence of such massive and compact 
star-forming galaxies at z = 3.2. 

4. DISCUSSIONS 
4.1. Selection Bias 

We note here on a possible selection bias introduced 
by our use of the [Oiii] emission line as an indicator 
of star-forming galaxies. When we use the [Oiil] line, 
the galaxy sample tend to be biased towards galax¬ 
ies with more extreme ISM conditions. It has been 
found that high-z star-forming ga laxies tend to have 
much higher excitation states (e.g. [Masters et al.l[20T^ 
Holden et al.ll20l5 iSteidel et al.l 12014 iShimakawa et al l 

20I4I) . iShimakawa et al.l (l2014f ) perform the NIR spec¬ 
troscopic observations of the Ha emitters at z=2.2 and 
2.5 associated to the two protocluster fields, and have 
shown that the [Oiii]/Ha ratios measured from the 
stacked spectra are ~1.0-3.0 in the stellar mass range of 
1O®-1O^^M0. The extreme ISM condition is expected to 
be a common feature among high-z star-forming galax¬ 
ies, and we expect that the [Oiii] emission line is an ap¬ 


propriate tracer of normal star-forming galaxies at high 
redshifts. 

The [Oiii] emitters may also be biased to less dusty 
galaxies compared to the Ha emitters, since the [Oiil] 
emission line (5007A) is located at the slightly shorter 
wavelength than the Ha emission line (6563A) and hence 
more strongly affected by dust extinction. However , 
adopting the extinction curve of ICalzetti et al.l (I 2000 D . 
the dust extinction at the wavelength of the [Oiii] line is 
only rsj 1.3 times larger than that at the wavelength of the 
Ha line. Considering that high-z star-forming galaxies 
tend to have high [Oiii]/Ha ratios as mentioned above, 
the effect of dust extinction for the [Oiii] line would not 
introduce a strong bias to less dusty galaxies. 

Moreover, metallicity of galaxies may also affect the 
strength of [Oiii] emission. Since lower metallicity 
leads to higher stellar temperature, [Oiii] line becomes 
stronger. Given the well known mass-metallicity rela¬ 
tion of star-forming galaxies (e.g. lErb et al]l2006aD . this 
metallicity effect may result in a possible bias towards 
lower stellar masses for [Oiii] emitters as compared to 
Ha emitters. 

In order to verify those selection bia ses, HiZELS (the 
High-redshift(Z) Emissio n Line Survey; iBest et aI1l201d : 
iSobral et al.l 12013112014D offers a very unique sample of 
dual emitters. They used a pair NB filters to capture 
[Oiii] and Ha emission lines at the same redshift, and 
constructed the samples of [Oiii] emitters and Ha emit¬ 
ters at z = 2.23. We will address the selection biases 
between the two samples based on this unique data sets 
in a forthcoming paper (D. Sobral, private communica¬ 
tion). 

4.2. Galaxy Growth from z = 3.2 to z = 2.2 

In Section 3.4, we show that there is no significant 
change in the location of the main sequence of star¬ 
forming galaxies between z = 3.2 (3.6) and z = 2.2 
(2.5), but the galaxy distributions on the sequence are 
different between the two epochs. In this section, we as¬ 
sume that the difference in galaxy distributions on the 
M*-SER plane between the [Oiii] emitters and Ha emit¬ 
ters is simply due to the evolution of star-forming galax¬ 
ies between the two epochs and discuss the stellar mass 
growth of galaxies from z = 3.2 to z = 2.2. From our re¬ 
sult that the location of the main sequence is unchanged 
during this time interval (IGyr), which is represented 
by SFR = as defined for the Ha emitters at 

z=2.2 and 2.5 (Figure [7]), we can put some constraints 
on the history of star formation and thus that of the stel¬ 
lar mass growth. In order to stay on the same main se¬ 
quence, the simplest evolutionary path would be that the 
individual star-forming galaxies evolve along the main 
sequence. This assumption should be valid if the galax¬ 
ies keep forming stars at the rates above our threshold 
of the Ha NB imaging, i.e. SFR > 4 [M 0 yr“^] (dust- 
uncorrected) and EWrest(Ha -I- [Nn]) > 40 [A]. 

The stellar mass growth between z = 3.2 and z = 2.2 
can be approximately tracked by the following derivative 
equation; 

dM^/dt = (I - i?) X SFR = (I - i?) X 129 (16) 

where the return mass fraction Ris ^ 0.3 for the Salpeter 
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IMF. Using this equation, a galaxy with M* = 10®Mq 
at z = 3.2 can increase their stellar mass by a factor of 
10 to ~ 1.1 X IO^^Mq by z = 2.2, while a galaxy with 
M* = 10 ^^Mq can grow in mass by a factor of 2. There¬ 
fore, more than 50% up to ^ 90% of the stellar mass of 
the star-forming galaxies at z = 2.2 can be formed dur¬ 
ing the 1 Gyr time interval between z = 3.2 and 2.2. The 
majority of the galaxies with M* > IO^Mq that we see at 
z = 3.2 would grow to massive galaxies of M* > IO^^Mq 
at z = 2.2, if they keep their high star formation activi¬ 
ties. 

Note also that, in this simple model, galaxies climb up 
the main sequence, i.e. SFR increases a lot from z = 3.2 
to z = 2.2 as the stellar mass grows. This indicates that 
the star formation activities of galaxies at z > 3 are ac¬ 
celerated towards the peak epoch of galaxy formation at 
z ^ 2. In this respect, z > 3 is the pre-peak epoch of 
galaxy formation. In order to achieve such an increasing 
star formation activity, an increasing rate of gas infall 
from outside is required, since otherwise the quick gas 
consumption would lower SFR as time progresses. In or¬ 
der to verify the presence of such continuous gas infall 
more quantitatively, we estimate the gas mass for the 
[Oiii] emitters from their SFR surface densities by assum - 
ing the Schmidt-Kennicutt relation (|Kenmcuttl 1 1 998bll . 
SFR surface densities (Esfr = SFR/TTr^) are estimated 
by using SFRs derived from UV luminosities in Section 
3.3 and the effective radius in Section 3.5. We then cal¬ 
culate the gas depletion time-scale of tdep = Mgas/SFR. 
The depletion time-scale of the [Oiii] emitters is mostly 
in the range of 0.2 — 0.4 Gyr, and shorter than 1 Gyr. 
This means that the [Oiil] emitters at z = 3.2 would 
consume all the remaining gas and terminate the star 
formation before z = 2.2 if there is no gas supply from 
the outside of galaxies. 

We have to mention that we have assumed the expo¬ 
nentially declining star formation history (SFH) in the 
form of SFR ^ exp(—t/r) in the SED fitting, while we 
now claim that the SFR increases with time from z = 3.2 
to z = 2.2. In order to verify the impact of assumed 
SFHs on the resulting physical quantities in the SED fit¬ 
ting, we re-estimate the stellar masses and SFRs of the 
[Oiii] emitters by assuming the exponentially increasing 
SFH. In the case of the increasing SFH, the estimated 
stellar masses vary by only factor of 0.9 — 1.3 for most 
of our sample, while the SFRs derived with Hy values 
from the SED fitting can increase by a factor of ~ 1.4. 
However, such a modest offset would be systematic and 
would apply to both the Ha and [Oiii] emitter samples. 
Therefore it should not change our results significantly. 

In reality, some galaxies would stop their star forma¬ 
tion and evolve to quiescent galaxies by z = 2.2, al¬ 
though this quenching process should happen on a rel¬ 
atively short timescale so that they do not significantly 
appear on the lower side of the main sequence and break 
its clear sequence. Also, we have ignored the effect of 
galaxy-galaxy mergers which can also increase the stellar 
mass of galaxies. Moreover, some galaxies would pop out 
all of a sudden on the main sequence with M* > IO^^Mq 
sometime between z = 3.2 and z = 2.2, which were below 
M* < IO^Mq at z = 3.2 or somewhere off the main se¬ 
quence. Those galaxies should form stars at even higher 
rates such as in a starburst mode, and the fraction of 


stars that are formed between the two epochs can be 
larger than 90%. 

The presence of those missing galaxies that are not 
considered in the simple model above is indicated by the 
comparison of number densities of the [Oiil] emitters at 
z = 3.2 and the Ha emitters at z = 2.2. The number 
density of the [Oiii] emitters at z=3.2 with M* > IO^Mq 
is 1.7 X 10“^ [Mpc“^], while that of the Ha emitters at 
z=2.2 with M* > lOi°M 0 is 2.7 x IQ-^ [Mpc'^j. Here 
we have taken into account the mass growth predicted 
by the above simple model. The latter number is ^1.6 
times larger. It suggests that, some galaxies may actually 
appear on the main sequence suddenly between z = 3.2 
and z = 2 . 2 , if we consider that there is no selection bias 
between the [Oiil] emitters and the Ha emitters (see 
Section 4.1). 

In any case, it is likely that star-forming galaxies grow 
at an accelerated pace during this time interval, assuring 
that this epoch is critically important for galaxy forma¬ 
tion. 

We also investigate the size growth of galaxies from z = 
3.2 to z = 2.2 by assuming that the mass-size relation 
is unchanged between z = 3.2 and z = 2.2 as suggested 
in Section 3.6. Using the mass-size relation of late -type 
galaxies at z ^ 2.75 from Ivan der Wei et al.l (I2014I1 . the 
effective radius of a galaxy with M^, = IO^Mq at z = 3.2 
would grow in size by a factor of ^ 1.5 by z = 2.2, and 
the size growth ratio does not depend much on the initial 
stellar mass of galaxies at z = 3.2. The size growth is not 
so strong from z=3.2 to z=2.2 as compared to the mass 
growth that we just discussed above. Considering the 
growth of the stellar mass and the size of galaxies from 
z = 3.2 to z = 2.2 together, we can also estimate the 
evolution in stellar mass surface density. It is predicted 
to grow by a factor of 5 for a galaxy with M^, = 10® Mq 
at z = 3.2. 

5. SUMMARY 

In this study, we construct an [Oiii] emitter sample at 
z > 3 in SXDF from the NB imag ing data taken wit h 
MOIRCS on the Subaru telescope (|Tadaki et al.ll201^ . 
We identify 27 and 7 [Oiii] emitters at z = 3.2 and 3.6, 
respectively. Some objects in our [Oiii] emitter sam¬ 
ple might be contributed by AGNs based on the rest- 
frame UVJ diagram and the Spitzer/MlPS detections. 
The spectroscopic observation is required to confirm the 
presence of AGNs and we do not exclude these objects 
in this study. Using the multi-wavelength data and HST 
high resolution images, we investigate their basic physi¬ 
cal properties, and compared them with those of the Ha 
emitters at z = 2.2 and 2.5 in the same field. 

• The stellar mass and the dust-extinction-corrected 
SFRuv of the [Oiil] emitters show a clear correla¬ 
tion as seen in other previous studies over a wide 
redshift range. Comparing our [Oiil] emitters at 
z = 3.2 and 3.6 with the H a emitters at z = 2. 2 
and 2.5 in the same field from iTadaki et ahl (j2013l) . 
the location of the M*-SFR relation of the [Oiil] 
emitters at z = 3.2 and 3.6 is almost the same as 
that of the Ha emitters at z = 2.2 and 2.5. 

• Although the location of the relation is almost 
the same between the [Oiii] and Ha emitters, the 
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galaxy distributions on the M*-SFR plane are dif¬ 
ferent in the sense that the [Oiil] emitters at z = 
3.2 and 3.6 tend to have lower stellar masses and 
SFRs as compared to the Ha emitters at z = 2.2 
and 2.5. 

• If we assume that the different galaxy distributions 
on the main sequence are due to the evolution of 
star-forming galaxies from z = 3.2 to z = 2.2, and 
that star-forming galaxies simply evolve along the 
constant star-forming main sequence in this time 
interval, galaxies with M* = 10®-10 ^^Mq can ob¬ 
tain ^ 90-50% of their stellar masses within just 
a Gyr from z = 3.2. Galaxies climb up the main 
sequence, and their star formation rates also in¬ 
crease a lot as their stellar masses grow. Although 
we consider only the simple model without outflows 
or mergers, we infer that galaxy formation activi¬ 
ties at z > 3 are accelerated towards its peak epoch 
at z 2. 

• We investigate the sizes of the [Oiii] emit- 
ters measured from th e HST H-band images 
(|van der Wei et al.ll2012tl . The size distribution of 
the [Oiil] emitters at z = 3.2 and 3.6 with respect 
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